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Group classification of a class of nonlinear fin equations is carried out exhaustively. Additional 
equivalence transformations and conditional equivalence groups are also found. These allow to 
simplify results of classification and further applications of them. The derived Lie symmetries 
, are used to construct exact solutions of truly nonlinear equations for the class under consider- 

ed ■ ation. Nonclassical symmetries of the fin equations are discussed. Adduced results amend and 

, essentially generalize recent works on the subject [M. Pakdemirli and A.Z. Sahin, Appl. Math. 

Lett., 2006, V.19, 378-384; A.H. Bokhari, A.H. Kara and F.D. Zaman, Appl. Math. Lett., 2006, 
V.19, 1356-1340]. 
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. 1 Introduction 

> ; 

Investigation of heat conductivity and diffusion processes leads to interesting mathematical 
models which can be often formulated in terms of partial differential equations. In general 
case coefficients of model equations explicitly include both dependent and independent model 

■ variables that makes difficulties in studying such models. 

. In this letter the class of nonlinear fin equations of the general form 

u t = (D(u)u x ) x + h(x)u, (1) 

| where D u ^ 0, is investigated within the symmetry framework. Here u is treated as the di- 

mensionless temperature, t and x the dimensionless time and space variables, D the thermal 
conductivity, h = —N 2 f{x), N the fin parameter and / the heat transfer coefficient. 

The condition D u = corresponds to the linear case of Q which was completely investigated 
with the Lie symmetry point of view long time ago pi HUj. Moreover, the sets of the linear and 
nonlinear equations of form Q can be separately investigated under restriction with point 
symmetries. That is why the linear case is excluded from consideration in the present letter. 

The problem of group classification for the degenerate case h = (i.e. the class of nonlinear 
one-dimensional diffusion equations) was first solved by Ovsiannikov [51 llUj. The equations of 
form Q with h being a constant are in the class of diffusion-reaction equations classified by 
Dorodnitsyn 013]. Group classification of the subclass where the thermal conductivity is a power 
function of the temperature was carried out in ^5j. We keep the above cases in presentation of 
results for reasons of classification usability. Note also that Lie symmetries of the class of quasi- 
linear parabolic equations in two independent variables, which has a wide equivalence group and 
covers all the mentioned classes, were classified in PJEJ- 

Recently Lie symmetries and reductions of equations from class (0) were considered in a 
number of papers f2\ IT2]. (See ibid for references on physical meaning and applications of 
equations Q.) In contrast to these papers, study in our letter is concentrated on rigorous and 
exhaustive group classification of the whole class Q and construction of exact solutions for truly 
nonlinear 'variable-coefficient' equations from this class. Additional equivalence transformations 
and conditional equivalence groups are also found. These allow to simplify results of classification 
and further applications of them. To find exact solutions, we apply both classical Lie reduction 
and nonclassical symmetry approaches. 
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2 Group classification and related problems 



Group classification of class (^Q) is performed in the framework of the classical approach ^U] . All 
necessary objects (the equivalence group, the kernel and all inequivalent extensions of maximal 
Lie invariance algebras) are found. Moreover, we extend the classical approach with additional 
equivalence transformations and conditional equivalence group for simplification of the classifi- 
cation results. 

The equivalence group G~ of class Q is formed by the transformations 



t = <5ii + <5 2 , 



6 3 x + 64, 



6 5 u, D = S^D, 



h = S^h, 



where 6i, i = 1,... ,5, are arbitrary constants, 616365 ^ 0. The connected component of the 
unity in G~ is formed by continuous transformations having 6± > 0, 63 > and £5 > 0. 
The complement discrete component of G~ is generated by three involutive transformations of 
alternating sign in the sets {t,D,h}, {x} and {u}. 

The kernel of the maximal Lie invariance algebras of equations from class (^Q) coincides with 
the one-dimensional algebra (dt). 

All possible C-inequivalent cases of extension of the maximal Lie invariance algebras are 
exhausted by ones adduced in Table 1, where 



h l (x) = eexp 



x 2 + p 



-dx 



; p € {—1,0, 1}, e = ±1, a 6 {0, 1} mod G~; n, q ^ 0. 



Table 1. Results of group classification 
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D(u) 


h(x) 


Basis of A max 
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dt 
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dt, d x 
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V 


x~ 2 


d t , 2td t + xd x 


4 


u n 


ex q 


d t , -qntdt + nxd x + (q + 2)ud u 


5 


u n 


ee x 


d t , -ntd t + nd x + ud u 


6 


u -4/3 


h}-(x) 


d t , -Aqtdt + 4(a; 2 + p)d x - 3(4x + q)ud u 


7 


V 





d t , d x , 2td t + xd x 


8 


(u+1)- 1 


£ 


d t , d x , e £t d t + £e £ \u + l)d u 


9 


e u 





d t , d x , 2td t + xd x , xd x + 2d u 


10 


u n , n ^ — | 


£ 


dt, d x , e~ ent (dt + £ud u ), nxdx + 2ud u 


11 


(u+a) n , n + -\ 





d t , d x , 2td t + xd x , nxd x + 2(u+a)d u 


12 


u -4/3 


£ 


dt, d x , ez et (dt + £ud u ), 2xd x — 3ud u , x 2 d x — 3xud u 


13 


(«+a)" 4 / 3 





dt, d x , 2tdt + xd x , 2xd x — 3(u+a)d u , x 2 d x — 3x(u+a)d u 



The Case 6 was missed in [TT] and in subsequent papers on the subject. The parameter- 
function h 1 equals to the following functions depending on values of p: 



P 



-1: h 1 



x - 1 



g/2 



X + l 

Additionally we can assume q 



p = 0: h 1 =£e~ q/x , p = l: h 1 
T mod G~ if p = 0. 



ee 



q arctan x 
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Some cases from Table 1 are equivalent with respect to point transformations which obviously 
do not belong to G~. These transformations are called additional equivalence transformations 
and lead to simplification of further application of group classification results. The pairs of 
point-equivalent extension cases and the corresponding additional equivalence transformations 
are 

6 p= o — * 5h=-4/3 '■ t = t, x = x , u = x u; 

x — 1 

6 p =_l -» ^n=-4/?,,q= q /2- * = t, X = — - ~, U = 2~ ?J ' 2 (x + lfu] 

ll a /o — * Ha=o, 13 Q ^o — > 13 a =o (« = — 3) ■ t = t, x = x, u = u + a; 
10 -» 110=0, 12 -» 13 Q=0 (n = -f ) : t = j^e £nt , x = x, u = e~ et u. 

The latter transformation was adduced e.g. in [3] . Case 6 with p = 1 is reduced to Case 4 only 
over the complex field. Note also that Case 8 is reduced by the similar transformation 

t = -\e~ £ \ x = x, u = e~ £ \u + 1) 

to the equation = (u^u^x — £ which is not in the class under consideration. 

There are no other additional equivalence transformations in class (|T|). Therefore, up to 
point equivalence, possible cases of extension of maximal Lie invariance algebras are exhausted 
by Cases 1-5, 6 p= i, 7-9, ll a =o and 13 a= o. 

The singularity of the diffusion coefficient D = u~ 4 / 3 with a number of different values of h 
admitting extensions of Lie invariance algebra can be explained in the framework of conditional 
equivalence groups. The equivalence group is extended under the condition D = ii~ 4//3 . More 
precisely, the equivalence group of the subclass of equations Q with D = -u -4 / 3 is formed 
by the transformations 

t = b]t + 02, x = - — , u = ±di(d5X + 06 J u, h = a% h, 

d 5 x + 6 

where Si, i = 1, . . . , 6, are arbitrary constants, 8\ > and ^3^6 — 6465 = ±1. G± is a non-trivial 
conditional equivalence group of class ((Hi. Two first additional equivalence transformations 
belong to G~. 

Another example of a conditional equivalence group in class Q arises under the condition 
h = 0. The equivalence group G~ of the whole class is then extended with translations with 
respect to u, i.e. the complete equivalence group of nonlinear diffusion equations (h = 0) is 
formed by the transformations 

t = 5\t + 62, x, = 5%x + 64, u = 5^u + 5q, D = 5i 1 5 2 D, 

where Si, i = 1,...,6, are arbitrary constants, S1S3S5 7^ 0. The third additional equivalence 
transformation belongs to G^- 

The subclass of equations (^Q) with h being a constant admits an extension of generalized 
equivalence group. The prefix "generalized" means that transformations of the variables t, x 
and u can depend on arbitrary elements jSJ. The associated generalized equivalence group G% 
is generated by transformations from G~ and the latter kind of additional equivalence transfor- 
mations, where e is replaced by h. 

Knowledge on conditional equivalence groups allows us to describe the set of admissible 
(from-preserving) transformations in class (^Q) completely. See e.g. ^H] and references therein. 

Note also that the subclass of equations (0) possessing nontrivial local conservation laws is 
exhausted by ones with h is a constant. Then the corresponding space CL of conservation laws 
is two-dimensional. The conserved vectors and the characteristics of basis elements of CL are 

(xe- ht u, e- ht (-xDu x +fD)), xe~ ht and (e _w «, -e~ ht Bu x ), e~ ht . 
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3 Similarity solutions 



Cases 7-13 of Table 1 are presented by 'constant coefficient' diffusion-reaction equations. More- 
over, all of these cases either are usual nonlinear diffusion equations or can be reduced to them 
by additional equivalence transformations. Exact solutions of 'constant coefficient' diffusion- 
reaction equations have been already investigated intensively. See for example, [3*1 14"[ |*"1 [Ti]. 
That is why we choose Cases 4-6 as representatives among truly nonlinear variable-coefficient 
fin equations, which are most interesting for Lie reduction. 
As shown in the previous section, the equation 

u t = {u- A ^u x ) x + h 1 {x)u (2) 

(Case 6 of Table 1) admits the two-dimensional (non-commutative) Lie invariance algebra q 
generated by the operators 

X 1 = d t , X 2 = -Aqtdt + 4(x 2 + p)d x - 3(4x + q)ud u . 

A complete list of inequivalent non-zero subalgebras of q is exhausted by the algebras (-X~i), (X 2 ) 
and(X 1 ,X 2 ). 

Lie reduction of equation Q to an algebraic equation can be made with the two-dimensional 
subalgebra {Xi,X 2 ) which coincides with the whole algebra g. The associated ansatz and the 
reduced algebraic equation have the form 

6.0. (X l5 X 2 ): u = C{x 2 +p)- 3 / 2 (h 1 {x)y 3/ \ C 4 / 3 = ^(q 2 + 16p). 

Substituting the solution C = ±^^(q 2 + 16p) 3//4 of the reduced algebraic equation into the 
ansatz, we construct the exact solution 

u = ±^-(q 2 + lQpf/\x 2 + p)- 3 /W(*)r 3/4 
8 

of equation Q. 

The ansatzes and reduced equations corresponding to the one-dimensional subalgebras from 
the optimal system are following: 

6.1. (Xi): u = (ip(uj))~ 3 , uj = x; 3<p ww = h l (uj)Lp~ 3 ] 

6.2. (X 2 ): u= ((x 2 + P ) 1 / 2 {h 1 {x)) l / A ^(uj)y 3 ,u J = th 1 (x); 

3q 2 u 2 ip U}U} + \q 2 u^pu - 3v?~ 4 ^ + -^(q 2 + 16p)(p - eip~ 3 = 0. 

The obtained reduced equations obviously have partial exact solutions which lead to the above 
exact solution of equation © and can be constructed via reduction to algebraic equations. The 
problem is to find some different solutions. We are only able to reduce the order of equation 
6.1. Namely, in the variables 

y = (u; 2 +p)- l l 2 {h l (u J ))- 1 /^, V> = (^ 2 +p)- l/2 {h\u)y l /\{u> 2 +p)<p u - u<p) 

constructed with the induced symmetry operator 4(u> 2 +p)d LU + (4o> + q)tpd,p equation 6.1 takes 
the form (4ijj — qy)ip y + qip + 4py = |ey~ 3 . A better way for construction of exact solutions for 
equations of Case 6 with p ^ is to map them to Cases 4 and 5 with additional equivalence 
transformations and then study the latter cases. 

Let us review results on Lie reduction of Cases 4 and 6. For each from these cases we 
denote the basis symmetry operators adduced in Table 1 by X\ and X 2 . Structure and list of 
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inequivalent subalgebras of the Lie invariance algebras are the same as ones in Case 6. The 
associated ansatzes and reduced equations have the form (e' = signi): 

4.0. (X 1 ,X 2 ): u = Cx 3 ¥, (q + 2)(nq + n + q + 2)C n+1 + en 2 C = 0; 

4.1. (Xi): u = (yj(o;))^T, u = x, (p uu + e(n + l)u^<p^ = ifn^-1; 

u = exp((p(u)), oj = x, ifuu + eujie'P = if n = —1; 

g+2 i .q + 2 .1 

4.2. (X 2 ): it = |i| "« u = \t Ux, (^ipuj)^ + ecu q <p + e' if - e ~uip u = 0; 

nq q 

5.0. (Xi,X 2 ): u = Cen, (n + l)C n+1 + en 2 C = 0; 

5.1. (Xi): u = (y)(o;))^+r, w = x, ip wti} +e(n + i)e u tp*+i = if n / -1; 

u = exp(</j(u;)), w = x, + ee^^ = if n = —1; 

5.2. (X 2 ): it = w = x + In |i|, (ip n ip u )u + ee^ip + e'n~ l (p - e'<p u = 0. 

Reduction to algebraic equations gives the following solutions of the initial equations: 

/ q + 2 . \ ~~ q+2 

4. u = 7r\ n< l + n + q + 2) \ x n ; 

\ en^ / 

/ n + 1\ ~™ * 

There are Emden-Fowler and Lane-Emden equations and their different modifications among 
the reduced ordinary differential equations. Solutions of these equations are known for a number 
of parameter values (see e.g. ^H])- As a result, classes of exact solutions can be constructed for 
fin equations corresponding to Cases 4 and 5 of Table 1 for a wide set of the parameters n and q. 



4 On nonclassical symmetries 

We also study conditional (nonclassical) symmetries of equations from class (0). As well-known, 
the operators with the vanishing coefficient of dt gives so-called 'no-go' case in study of condi- 
tional symmetries of an arbitrary (1 + l)-dimensional evolution equation since the problem on 
their finding is reduced to a single equation which is equivalent to the initial one (see e.g. [TB]). 
Since the determining equation has more independent variables and, therefore, more freedom 
degrees, it is more convenient often to guess a simple solution or a simple ansatz for the deter- 
mining equation, which can give a parametric set of complicated solutions of the initial equation. 
For example, the fin equation 

u t = (u~ l u x ) x + xu (3) 

is conditionally invariant with respect to the operator d x + tud u . The associated ansatz u = 
e tx (p(u>), u = t, reduces equation © to the equation ip u = 0, i.e. u = Ce tx is its non-Lie exact 
solution which can be additionally extended with symmetry transformations. 

It is known that non-linear diffusion equations (i.e. equations (^Q) with h = 0) possess condi- 
tional symmetry operators which have non-vanishing coefficients of dt and are inequivalent to 
Lie invariance operators only in case of exponential diffusion coefficients (Case 9 of Table 1 and 
equivalent equations). Moreover, solutions associated with such conditional symmetry operators 
are still Lie invariant. The opposite situation is for h 7^ 0. There exist conditional symmetry 
operators of equations which have non- vanishing coefficients of dt, are inequivalent to Lie 
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invariance operators and even lead to truly non-Lie exact solutions. Thus, consider again equa- 
tion (0). It admits also the conditional symmetry operator dt + xud u . The associated ansatz 
u = e tx ip(uj), uj = x, reduces equation © to the equation (tp^ 1 ^^)^ = 0. The general solution 
ip = C\e C2X of the reduced equations gives a solution of (jHJ), which is simplified to the above 
constructed one with symmetry transformations. 

Exhaustive description of nonclassical symmetry operators of equations (0) will be a subject 
of a forthcoming paper. 
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